Confinement and the quark Fermi-surface in SU(2N) QCD-like theories 
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Yang-Mills theories with a gauge group SU{Nc 7^ 3) and quark matter in the ftindamental rep- 
resentation share many properties with the theory of strong interactions, QCD with Nc = 3. We 
show that, for Nc even and in the confinement phase, the gluonic average of the quark determinant 
is independent of the boundary conditions, periodic or anti-periodic ones. We then argue that a 
Fermi sphere of quarks can only exist under extreme conditions when the centre symmetry is spon- 
taneously broken and colour is liberated. Our findings are supported by lattice gauge simulations 
for Nc = 2 ... 5 and are illustrated by means of a simple quark model. 
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In particle physics, understanding the theory of strong 
interactions, i.e., QCD, under extreme conditions is key 
for a successful description of exotic matter which fea- 
tures in the evolution of the universe or cold compact 
stars. The property of QCD which is most relevant for 
shaping the QCD phase diagram (if presented as a func- 
tion of temperature and matter density) is colour con- 
finement. Under normal conditions, confinement implies 
that quark matter is organised in terms of hadrons only, 
and quarks and gluons are only part of the particle spec- 
trum under extreme conditions. 

By means of lattice gauge simulations, a precise picture 
of hot QCD matter at small densities has emerged over 
the last two decades: Centre symmetry is spontaneously 
broken in the hot de-confinement phase for temperatures 
above a certain critical value. Colour is liberated and, 
in a theory with heavy quarks only, the Polyakov line 
expectation value serves as an order parameter [l|, . 

On the other hand, the situation at high densities and 
small temperatures is far from being clear. The reason is 
the lack of first principle QCD results which would help 
to scrutinise the proposals for the properties of matter in 
this regime. Lattice gauge simulation techniques cannot 
be applied because of the severity of the so-called sign 
or overlap problem. These problems are absent for the 
so-called two-colour QCD, and intriguing results, even 
for large quark chemical potential, have been accumu- 
lated over the recent years For an investigation 
of cold and dense matter in SU{Nc > 3) QCD(-like) 
theories, we are still awaiting major conceptual achieve- 
ments. Promising recent attempts abandon standard lat- 
tice Monte-Carlo techniques and are based upon stochas- 
tic quantisation [6] or worldline numerics [7|. 

Perturbative QCD and QCD inspired quark models 
have been a valuable tool for revealing mechanisms which 
might operate in the cold and dense phase of QCD. On 
this basis, the QCD phase diagram has gained a lot of 
renewed interest when findings suggested that its struc- 
ture is far more complex than it had been suggested 
for decades: at the highest densities, it is expected that 



quark matter is forming a colour-superconductor which 
carries along a rich phase structure on its own (for a re- 
cent review see Q). Studies of the Gross-Neveu model 
indicate that dense QCD (at low temperatures) might 
form an inhomogeneous baryonic crystal [9l-[Tl|. 

Employing arguments based upon the large Nc expan- 
sion, it has been recently suggested that the tight rela- 
tion between confinement and spontaneous chiral sym- 
metry breaking (inherent for zero density QCD) gets al- 
leviated [13, [iJl ■ In the phase diagram as a function of 
chemical potential and temperature, the phase boundary 
for chiral restoration might deviate from the boundary for 
deconfinement. Most interesting, a phase for which con- 
finement is still intact while chiral symmetry is restored 
has attracted a lot of interest. This so-called quarkyonic 
phase is characterised by a Fermi sphere of quarks while 
the outer shell of the Fermi sphere necessarily consists of 
baryons because of confinement [T^ . 

In this paper, we point out that the properties of dense 
quark matter in SU (Nc) QCD-like theories are vastly dif- 
ferent depending on whether the number of colours, Nc, 
is even or odd. In the confining phase, the quark deter- 
minant is averaged over centre-transformed gluonic back- 
ground fields. Underpinned by lattice gauge simulations, 
we will show that this averaged quark determinant is in- 
sensitive to the boundary conditions of the quarks, peri- 
odic or anti-periodic ones. Since the formation of a quark 
Fermi sphere is crucially linked to anti-periodic bound- 
ary conditions (see below for an illustration by a quark 
model), we will argue that a Fermi sphere of quarks can 
only exist in the deconfined phase of SU{Nc) QCD-like 
theories with Nc being even. 

QCD-like theories are SU{Nc) Yang- Mills theories cou- 
pled to fermions ( "quarks" ) in the fundamental represen- 
tation. We here adopt the lattice regularisation based 
upon a toroidal space-time lattice with lattice spacing a 
and extension Nt x N^. The gluonic degrees of freedom 
Uf_i{x) e SU{Nc) satisfy periodic boundary conditions, 
e.g.,Uf_i{xo+Nta,x) = C/^(a;o, x) . Quark fields (/(a;) are 
associated with the lattice sites. Because of the Fermi 



2 



statistics of the bare quark fields, these fields satisfy an- 
tiperiodic boundary conditions: 

q{xo + Nta,x) = (-l)g(xo,f). 

Using the Wilson action for the gluonic fields, the parti- 
tion function is given by 

Z = y"2?J7^ DetAAf[t/] e^^-.f-- ^^-f^), (1) 

where /3 = 2Nc/g'^ is given in terms of the Yang-Mills 
gauge coupling g, and P^^i, is the standard plaquette. 
The determinant in ([T]) arises from the integration over 
the quark fields. The subscript "A" indicates that the 
quarks were subjected to antiperiodic boundary condi- 
tions. Here, we work with Wilson quarks where 



M[U] 



^ 4 

(m + A)6^y - ^^[{l- 7m) Ut,{x) (5:r+M,i 



(2) 



where 7^ are the Hermitian Dirac matrices and m is the 
current quark mass in units of the lattice spacing. For 
an investigation of the quark determinant, we insert 



1 = J dQ s(q ^BetAMp] 
into ([1]) to write the partition function as 



2 ^ J dQQPAiQ), (3) 

Pa{Q) = JvU^s(^Q-BetAM[U]^e^^"^^\ (4) 

where Pa{Q) is the probability distribution of the quark 
determinant. Note that Pa{Q) can be calculated using 
Monte-Carlo techniques for pure Yang-Mills theory. Note 
that the probability distribution of the quark determi- 
nant of full QCD, Pfuii(Q), is related to that of pure 
Yang-Mills theory, i.e., Pa{Q), by Pfuii(Q) = QPa{Q). 
In particular for large lattices, little statistics is expected 
for the large Q regime which is more relevant for the 
simulation of Yang-Mills theory with dynamical quarks 
included. If this regime is under considerations, it is ad- 
visable to include the determinant in the simulation 
and to use refined simulation techniques such as those 
m , [1^ when finite temperatures and densities are ad- 
dressed. The intermediate Q regime will turn out to be 
sufficient to illustrate our findings below, and thus only 
quenched simulations are used throughout this paper. 

Let us now consider a Roberge- Weiss transforma- 
tion [17[ in the gluonic functional integral in For 
a fixed xq = t, we consider 

Uo(t,x) ZnUo{t,x) Vx, (5) 
r27ri 



Given the invariance of the gluonic action and of the Haar 
measure VU^^ we find: 

Pa{Q) = y"p[/^<5(Q-Det^M[z„C/])e^™[^l. (7) 

Reintroducing the quark fields for a moment, 

DetAM[z„C/]= fvqVqcxp{j2q{x)M^yq{y)y (8) 

xy 

we explore the virtue of the transformation: 

U^ix) = Qix)U^{x)QHx + fi), (9) 
q^'ix) = nix)qix) , (10) 
Q{xo,x) — z„ for xo > t , r2(a;o, x) = 1 else. (11) 

Note that t has been defined above ^ . We point out that 
fl{x) satisfies the boundary condition D,{xo + Nta,x) = 
Zn^l{xQ,x). The transformation ([9]) does not change the 
boundary conditions of the link fields and leaves the 
Yang-Mills action invariant. It can be therefore consid- 
ered as a gauge transformation of pure SU{Nc) Yang- 
Mills theory. Note, however, that the transformation ([TU]) 
does change the boundary conditions of the quark fields: 

q'^{xo + Nta,x) = (-1) Zn q" {xq , x) . (12) 

It therefore does not qualify as a gauge transformation of 
the full theory. It, however, reveals an important prop- 
erty of the distribution Pa{Q)- Specialising to an even 
number of colours and choosing n = iVc/2, z„ = —1, we 
obtain from (gl) using PfTTj) : 



BetAM[znU] = BetpM[U] 



(13) 



exp 



{^n}, 0<n<N,-l. (6) 



where the subscript "P" signals that the quark fields now 
obey periodic boundary conditions. Inserting the last 
result (IT3I) into ([7]) and using the gauge invariance of 
Yang-Mills action and Haar measure, we finally obtain: 

Pa{Q) = J VUf,6(Q-DctpM[U]^e^™^^^ = Pp{Q). 

In QCD-like theories for an even number of colours and, 
at least, for a finite volume (see comment below), the 
partition function is independent of the choice of the 
quark boundary conditions, periodic or antiperiodic ones. 
In the infinite volume limit, the centre symmetry corre- 
sponding to the transformation @ is not always realised: 
it has been known for a long time [H [3] that this sym- 
metry is spontaneously broken under extreme conditions, 
high temperature and/or fermion densities. While peri- 
odic boundary conditions are associated with an insta- 
bility of the partition function for vanishing temperature 
corresponding to Bose-Einstein condensation, antiperi- 
odic boundary conditions are the essential ingredient for 
building up a Fermi sphere in the dense phase. Hence, we 
argue that in the confining phase of SU{2N) QCD-like 
theories, a quark Fermi surface is unlikely to exist. We 
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FIG. 1: Probability distribution for Inr in comparison for SU(2) and SU(3) (left) and for SU(4) and SU(5) (right). 



stress, however, that the deconfinement phase at high 
densities (and low temperatures) might weh feature a 
Fermi sphere of quarks. 

We are now going to iUustrate our findings using lattice 
gauge simulations as well as a simple quark model. We 
have carried out simulations for SU(A^c), Nc = 2,3,4,5 
using a 4^ space time lattice and a Dirac mass m = 0.01. 
The Wilson parameters /3 have been chosen such that in 
all cases the lattice spacing a in units of string tension 
a was roughly constant, i.e., cra^ = 0.467(10). Keeping 
aa^ fixed when the number of colours Nc is increased 
also complies with the so-called 't Hooft limit where 
g'^Nc =constant. The simulation parameters are sum- 
marised in table HI 



group 
1^ ■ N 

-/n ai VP ♦ C 



SU(2) 
2.2 
3.63 



SU(3) 
5.63 
3.2 



SU(4) SU(5) 
10.99 16.57 
2.9 3.0 



TABLE I: Simulation parameters; 9902 configurations were 
used for SU(2) and SU(3), 8000 for SU(4) and 4500 for SU(5) 
to estimate the expectation values. 

The determinants have been calculated exactly using 
the standard LU-decomposition. In order to explore the 
sensitivity of the quark determinants to the boundary 
condition, we define 



V)etpM[U] 
V)etAM[U = 1] 



y = 



Det aM[C/ = 1]' 



(14) 



and consider the ratio r[U] := x/y for a given lattice 
configuration. For even, we find with the results above 
that 



r [z„{/] - l/r[U] 



(r) = (1/r) 



(15) 



if the centre-symmetry is realised. Our numerical find- 
ings for these expectation values are summarised in ta- 
ble HI 



SU(2) SU(3) SU(4) SU(5) 

(r) 1.66 ± 0.02 1.365 ± 0.009 2.9977 ± 0.07 1.0469 ± 0.005 
(l/r-> 1.65 ± 0.02 1.587 ±0.02 2.87 ± 0.07 1.0571 ± 0.007 

TABLE II: Sensitivity parameters (r) and (1/r) for several 
gauge groups. 



We here find a clear coincidence between (r) and 
for the gauge group SU(2) while the corresponding pa- 
rameters for SU(3) are significantly different. We also 
observe a tendency that the difference fades away for 
increasing number of colours Nc- We finally present 
the probability distribution for the variable Inr in fig- 
ure [T] For gauge groups with Nc even, we expect that 
the corresponding histograms are symmetric with respect 
to In r — > — In r (see 1151) . This expectation is nicely con- 
firmed for the gauge groups SU(2) and SU(4). A clear 
asymmetry is observed for SU(3) while the asymmetry is 
very small for SU(5) if present at all. We also empirically 
observe that the width of the probability distribution de- 
creases for an increasing number of colours. 

Let us finally illustrate by means of a quark model 
how centre sector tunneling eliminates the quark Fermi 
surface. The theory which we are proposing is essentially 
a free theory of quarks which, however, interact with a 
constant centre background field. We here work in the ab 
initio continuum formulation. Parameterizing the centre 
background field by [H = diag(l, ...,1,1 — Nc)/Nc from 
the Cartan algebra) 



(16) 



(where T is the temperature) the partition function of 
our model is in Euclidean space 



J VqVq exp{q{i^ + {A„ + iii)"fo + im)q] 



J VqVq exp{g(j^ + im)q} 



(17) 
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where m is the quark mass and /i is quark chemical poten- 
tial. Thereby, p„ is the probability that the centre sector 
n is attained. In the quenched approximation, all cen- 
tre sectors occur with equal probability, i.e., Pn = ^/Nc- 
Note, however, that in a more QCD relevant setting dy- 
namical quarks induce a bias towards the trivial centre 
sector, i.e., po > Pn^o- Furthermore, An acts as constant 
gauge field which can be eliminated from the action at 
the expense of changing the boundary conditions for the 
quark fields. The partition function can be calculated in 
closed form: 

^ - E?'" n (l + -.e-[^(^^-^l/T) + ^te-[i^(.l+.l/T) ^ 

n p^q 

where z„ is related to the trace of the Polyakov line P 
line by 

P = exp{iA„/r}, trP = z„ 

Moreover, the product extends over spatial momenta p 
and E{p) = {m?+p)^/'^ . Expanding the brackets, because 
of the sum over the centre elements, only states with van- 
ishing A'^-ality contribute to the partition function assert- 
ing confinement. We here focus on the BEC instability. 
We consider temperatures T which are small compared 
to the mass gap implying that we may neglect the con- 
tribution of antiquarks to the free energy for /i > m. The 
free energy is then given by In Z w In p„ /?„ with 

Pn — exp|-^ j dE E \/ E"^ - im? In ^1 + z„e '^)}) 

where V is the spatial volume. Because of the spin of the 
quarks, we have used two states per Matsubara mode. 
For the search for Fermi surface effects, it is most in- 
structive to study the baryon number density: 

T dlnZ 1 /"^ 

b =- —— = — dEE Ve^^ p{E,T,fi), 

V dp, TT^ 7„ 

where p is defined by 

n ^ 

with the weights Wn = Pn Pn/ J2iPi Pi- decon- 
fined phase, tunneling between centre sectors stops be- 
cause of spontaneous symmetry breaking (on top of the 
explicit breaking by the quark determinant). The sum 
over the centre elements collapses to the trivial centre el- 
ement Po = 1, Pn^o = implying wq ^ 1, Wn^o = 0. 



In this case, p{E, T, p) can be interpreted as spectral 
density. This is given by the familiar Fermi function 
PdcconiE, T,p) = [e^E-t^y^ + , which features a 

Fermi surface for E k, p. 

Let us now consider the confinement phase of a SU(2) 
gauge theory. Despite of the bias towards the trivial cen- 
tre sector, we have po < \ because of centre sector tun- 
neling. The crucial observation is that while increasing 
the chemical potential p to approach the mass m from 
below, pi develops a singularity while po is perfectly fi- 
nite. This implies that the weights are given by wq = 0, 
wi = 1 for p — > m. Thus, p{E, T, p) is approximately 
given by 

and, hence, signals the BEC instability for p ^ m. 
Note that the sign of p is dictated by the centre element 
zi = —1. Because of this sign, the contribution to the 
baryon number density is negative and arises from centre 
dressed quarks. This contribution is genuinely different 
from the contribution from bare anti-quarks which can 
be neglected for the present choice of parameters. 

In conclusions, we have shown that the gluonic aver- 
age of the quark determinant of SU {2N) gauge theories 
does not depend on the type of boundary conditions for 
the quark fields as long as the centre symmetry is re- 
alised. Lattice gauge simulations for the gauge groups 
SU{2 ... 5) corroborate these findings. Our results sup- 
plement the recent findings of [ISi] where SU(2) deconfine- 
ment was brought in line with the sensitivity of the quark 
determinant to the boundary conditions. Our results 
may have far reaching phenomenological implications: (i) 
they exclude the (pre)formation of a quark Fermi sphere 
at finite densities in the confinement phase; (ii) since 
quarks are dressed with centre fields of the gluonic back- 
ground they escape the spin-statistics connection. This 
is in analogy to the anyons in solid state physics [Toj . 
While anyons only occur in 2+1 dimensions, our model 
is the first one of its kind which evades the spin-statistics 
connection in 3-1-1 dimensions, (iii) Exotic states of mat- 
ter, such as a Bose-Einstein condensate of quarks, might 
exist prior to deconfinement induced by density. 
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